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Abstract. Let A* > denote the largest possible value of A such that 
A 2 u = Xe u in B 

du 

u = — = on dB 

dn 

has a solution, where B is the unit ball in M. N and n is the exterior unit normal 
vector. We show that for A = A* this problem possesses a unique weak solution 
u* . We prove that u* is smooth if N < 12 and singular when N > 13, in which 
case u*(r) = -41ogr + log(8(A r - 2)(N - 4)/A*) + o(l) as r -> 0. We also 
consider the problem with general constant Dirichlct boundary conditions. 



1. Introduction 
We study the fourth order problem 



(1) 



A 2 u = Ae" in B 
u = a on dB 

^ = b on dB 



dn 

where a, b e R, B is the unit ball in R N , N > 1, n is the exterior unit normal 
vector and A > is a parameter. 

Recently higher order equations have attracted the interest of many researchers. 
In particular fourth order equations with an exponential non-linearity have been 
studied in 4 dimensions, in a setting analogous to Liouville's equation, in [3, 12, 24] 
and in higher dimensions by [1, 2, 4, 5, 13]. 

We shall pay special attention to (1) in the case a = b = 0, as it is the natural 
fourth order analogue of the classical Gelfand problem 

Au = \e u in Q 
u = on dfl 



(2) 



(0 is a smooth bounded domain in Mr) for which a vast literature exists [7, 8, 9, 
10, 18, 19, 20, 21]. 

From the technical point of view, one of the basic tools in the analysis of (2) is 
the maximum principle. As pointed out in [2], in general domains the maximum 
principle for A 2 with Dirichlct boundary condition is not valid anymore. One of the 
reasons to study (1) in a ball is that a maximum principle holds in this situation, 
see [6] . In this simpler setting, though there are some similarities between the two 
problems, several tools that are well suited for (2) no longer seem to work for (1). 
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As a start, let us introduce the class of weak solutions we shall be working with : 
we say that u G H 2 (B) is a weak solution to (1) if e" G L 1 (_B), u = a on dB, 



§^ = bondB and 

on 

[ AuAtp = \( e>, for all ip G C™(B). 
Jb Jb 

The following basic result is a straightforward adaptation of Theorem 3 in [2] . 

Theorem 1.1. ([2]) There exists A* such that ifO<\<\* then (1) has a minimal 
smooth solution u\ and if X > A* then (1) has no weak solution. 

The limit u* = lim>,/'A* u\ exists pointwise, belongs to H 2 (B) and is a weak 
solution to (1). It is called the extremal solution. 

The functions u\, < A < A* and u* are radially symmetric and radially de- 
creasing. 

The branch of minimal solutions of (1) has an important property, namely u\ is 
stable in the sense that 



(3) / (A^>A / e"V, V<p€C%>(B) t 

J B J B 

see [2, Proposition 37]. 

The authors in [2] pose several questions, some of which we address in this work. 
First we show that the extremal solution u* is the unique solution to (1) in the class 
of weak solutions. Actually the statement is stronger, asserting that for A = A* 
there are no strict super-solutions. 

Theorem 1.2. // 

(4) v G H 2 (B), e v G L\B), v\ 9B = a, §^\ dB < b 
and 

(5) f AvA<p>\* [ e"ip V^GC °°(B), </? > 0, 
Jb Jb 

then v — u* . In particular for A = A* problem (1) has a unique weak solution. 

This result is analogous to work of Martel [19] for more general versions of (2) 
where the exponential function is replaced by a positive, increasing, convex and 
superlinear function. 

Next, we discuss the regularity of the extremal solution u*. In dimensions = 
5, . . . , 16 the authors of [2] find, with a computer assisted proof, a radial singular 
solution U a to (1) with a = b = associated to a parameter A CT > 8(N — 2) (AT — 4). 
They show that X a < A* if N < 10 and claim to have numerical evidence that 
this holds for < 12. They leave open the question of whether u* is singular in 
dimension N < 12. We prove 

Theorem 1.3. If N < 12 then the extremal solution u* of (1) is smooth. 

The method introduced in [10, 20] to prove the boundedness of u* in low di- 
mensions for (2) seems not useful for (1), thus requiring a new strategy. A first 
indication that the borderline dimension for the boundedness of u* is 12 is Rellich's 
inequality [23] , which states that if N > 5 then 



N 2 (N-4) 2 f ip 2 

16 VR 4 



(6) (Ay) 2 > \' frr v^c »r), 



STABLE SOLUTIONS FOR A 2 u = Xe 



3 



where the constant TV 2 (TV — 4) 2 /16 is known to be optimal. The proof of Theo- 
rem 1.3 is based on the observation that if u* is singular then X*e u ~ 8(TV — 2)(TV— 
4)|a;|- 4 near the origin. But 8(TV - 2)(TV - 4) > TV 2 (TV - 4) 2 /16 if N < 12 which 
would contradict the stability condition (3). 

In view of Theorem 1.3, it is natural to ask whether u* is singular in dimension 
TV > 13. If a = b — 0, we prove 

Theorem 1.4. Let TV > 13 and a = b = 0. Then the extremal solution u* to (1) 
is unbounded. 

For general boundary values, it seems more difficult to determine the dimensions 
for which the extremal solution is singular. We observe first that given any a, b G R, 
u* is the extremal solution of (1) if and only if u* — a is the extremal solution of 
the same equation with boundary condition u = on dB. In particular, if A* (a, b) 
denotes the extremal parameter for problem (1), one has that A* (a, b) — e~ a A*(0, b). 
So the value of a is irrelevant. But one may ask if Theorem 1.4 still holds for any 
TV > 13 and any b G R. The situation turns out to be somewhat more complicated : 

Proposition 1.5. 

a) Fix TV > 13 and take any a G R. Assume b > —4. There exists a critical 
parameter b max > 0, depending only on TV, such that the extremal solution 
u* is singular if and only if b < b max . 

b) Fix b > —4 and take any a G R. There exists a critical dimension N mm > 
13, depending only on b, such that the extremal solution u* to (1) is singular 
if N > N mm . 

Remark 1.6. 

• We have not investigated the case b < —4. 

• If follows from item a) that for b G [—4, 0], the extremal solution is singular 
if and only if N > 13. 

• It also follows from item a) that there exist values ofb for which N mm > 13. 
We do not know whether u* remains bounded for 13 < N < jV mm . 

Our proof of Theorem 1.4 is related to an idea that Brezis and Vazquez ap- 
plied for the Gclfand problem and is based on a characterization of singular energy 
solutions through linearized stability (see Theorem 3.1 in [8]). In our context we 
show 

Proposition 1.7. Assume that u G H 2 (B) is an unbounded weak solution of (1) 
satisfying the stability condition 

(7) A f eV < / (A^) 2 , G C °°(S). 

J B J B 

Then A = A* and u = u* . 

We do not use Proposition 1.7 directly but some variants of it - see Lemma 2.6 
and Remark 2.7 in Section 2 - because we do not have at our disposal an explicit 
solution to the equation (1). Instead, we show that it is enough to find a sufficiently 
good approximation to w*. When TV > 32 we are able to construct such an ap- 
proximation by hand. However, for 13 < N < 31 we resort to a computer assisted 
generation and verification. 
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Only in very few situations one may take advantage of Proposition 1.7 directly. 
For instance for problem (1) with a = and b = —4 we have an explicit solution 



associated to A = 8(N — 2)(N — 4). Thanks to Rellich's inequality (6) the solution 
u satisfies condition (7) when N > 13. Therefore, by Theorem 1.3 and a direct 
application of Proposition 1.7 we obtain Theorem 1.4 in the case b = —4. 

In [2] the authors say that a radial weak solution u to (1) is weakly singular if 



For example, the singular solutions U a of [2] verify this condition. 
As a corollary of Theorem 1.2 we show 

Proposition 1.8. The extremal solution u* to (1) with b > —4 is always weakly 
singular. 

A weakly singular solution either is smooth or exhibits a log-type singularity at 
the origin. More precisely, if u is a non-smooth weakly singular solution of (1) with 
parameter A then (see [2]) 



In Section 2 we describe the comparison principles we use later on. Section 3 is 
devoted to the proof of the uniqueness of u* and Propositions 1.7 and 1.8. We prove 
Theorem 1.3, the boundedncss of u* in low dimensions, in Section 4. The argument 
for Theorem 1.4 is contained in Section 5 for the case TV > 32 and Section 6 for 
13 < N < 31. In Section 7 we give the proof of Proposition 1.5. 

Notation. 

• Br: ball of radius R in M. N centered at the origin. B = B\. 

• n: exterior unit normal vector to Br 

• All inequalities or equalities for functions in L p spaces are understood to 



u(x) 



41og|x| 



lim ru'(r) exists. 



lim u(r) +4 log r = log 



8(iV-2)(iV-4) 
A 



lim ru' (r) = —4. 



be a.e. 



2. Comparison principles 



Lemma 2.1. (Boggio's principle, [6]) If u £ C a {Br) satisfies 

( A 2 u > in B R 




u = — — =0 on 8Br 



then u > in Br. 
Lemma 2.2. Let u £ L 1 



(Br) and suppose that 




for all tp £ C 4 (Br) such that ip > in Br, <p\dB 
Br. Moreover u = or u > a.e. in Br. 



R 



= Then u > in 



For a proof see Lemma 17 in [2]. 
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Lemma 2.3. If u G H 2 (Br) is radial, A 2 u > in Br in the weak sense, that is 

/ AuA^>0 V(peC%°(B a ),<p>0 
Jb r 

and u\o Br > 0, §^\dB R < then u > in B R . 

Proof. We only deal with the case R — 1 for simplicity. Solve 

A 2 wi = A 2 u in Bi 

dui o «D 
ui = — — = on oBi 

on 

in the sense u\ G Hq(Bi) and f g AuiAip = J B AuAip for all ip G C§°(Bi). Then 
ui > in B\ by Lemma 2.2. 

Let U2 = u — ui so that A 2 it2 = in B\. Define / = Ai*2. Then Af = in B\ 
and since / is radial we find that / is constant. It follows that U2 — ar 2 + b. Using 
the boundary conditions we deduce a + b > and a < 0, which imply U2 > 0. □ 

Similarly we have 

Lemma 2.4. If u E H 2 (Br) and A 2 u > in Br in the weak sense, that is 
[ AuAif > Vifl G C$°(B R ), ip>0 

J B R 

and u\qb r = 0, §^\ob r < t/ien u > m Sr. 

The next lemma is a consequence of a decomposition lemma of Moreau [22] . For 
a proof sec [14, 15]. 

Lemma 2.5. let u G Hq(Br). Then there exist unique w,v G Hq(Br) such that 
u = w + v, w > 0, A 2 v < in Br and J Br AwAv = 0. 

We need the following comparison principle. 

Lemma 2.6. let u\, 1*2 G H 2 (B R ) with e Ul , e U2 G ^(Br). Assume that 

A 2 m < \e Ul in B R 

in the sense 

(8) / A Ul Aip<\{ e Ul ip V(peCZ°(B R ),ip>0, 

J B R J B R 

and A 2 U2 > Ae" 2 in Br in the similar weak sense. Suppose also 

dui du 2 , 

u\\dB R = u 2 \dB R and -^—\dB R = ^—\dB R - 

on on 

Assume furthermore that u\ is stable in the sense that 

(9) a / e"V</ (Aip) 2 , V^eC^(Br). 



Then 

ui < U2 in Br. 
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Proof.Let u = Ui — u%. By Lemma 2.5 there exist w,v £ Hq(Br) such that 
u = w + v, w > and A 2 v < 0. Observe that v < so w > u\ — u 2 . 
By hypothesis we have for all <p £ C^°(B R ), p> 0, 



A(ui - u 2 )Atp < A / (e ui - e U2 )p < A / (e ui - e U2 )p 

Br J Br J B R n[ui>U2] 

and by density this holds also for w: 
(10) / (Aw) 2 = / A( Ul -u 2 )Aw 

JBr JBr 



< A / (e" 1 - e U2 )w = A / (e" 1 - e U2 )w, 

J B R n[u 1 >u 2 ] JBr 

where the first equality holds because J B AwAv = 0. By density we deduce from 
(9): 

(11) A f e u 'w 2 < f (Aw) 2 . 

JBr JBr 

Combining (10) and (11) we obtain 

e Ul w 2 < f [e Ul -e U2 )w. 

Br JBr 

Since u\ — u 2 < w the previous inequality implies 

(12) < / (e Ul - e U2 - e Ul ( Ul - u 2 ))w. 

JBr 

But by convexity of the exponential function e" 1 — e U2 — e Ul (ui ~ u 2 ) < and we 
deduce from (12) that (e Ul — e U2 — e" 1 (ui — u 2 ))w = 0. Recalling that u\ — u 2 < w 
we deduce that u\ < u 2 . □ 

Remark 2.7. The following variant of Lemma 2.6 also holds: 

Let ui, u 2 £ H 2 (B R ) be radial with e ui , e" 2 G L 1 (B R ). Assume A 2 u x < Ae" 1 
in B R in the sense of (8) and A 2 u 2 > Xe U2 in Br. Suppose ui\gB R < u 2 \qb r and 
T^HsBk > ^^\dB R cmd that the stability condition (9) holds. Then u\ < u 2 in Br. 

Proof. We solve for u e H§(B R ) such that 

/ AuAip = I A( Ul -u 2 )Aip V^6(7 M (Bs). 

JBr JBr 

By Lemma 2.3 it follows that u > u\ ~ u 2 . Next wc apply the decomposition of 
Lemma 2.5 to u, that is u = w + v with w, v £ Hq(B r ), w>0, A 2 v < in B R and 
J Br AwAv = 0. Then the argument follows that of Lemma 2.6. □ 

Finally, in several places we will need the method of sub and supersolutions in 
the context of weak solutions. 

Lemma 2.8. Let A > and assume that there exists u £ H 2 (B R ) such that 
e^eL^Bn), 

I AuAp > A / e> for all p £ C^{B R ) 7 p > 
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9u , ^„ 
u = a, — < b on oB\. 
on 

Then there exists a weak solution to (1) such that u < u. 
The proof is similar to that of Lemma 19 in [2]. 

3. Uniqueness of the extremal solution: proof of Theorem 1.2 

Proof of Theorem 1.2. Suppose that v £ H 2 (B) satisfies (4), (5) and v ^ u* 
Notice that we do not need v to be radial. 
The idea of the proof is as follows : 



Step 1. The function 



u o = t:{u*+ v) 



2 

is a super-solution to the following problem 

A 2 u = X*e u + (ir)e u in B 



(13) 



u = a 



on dB 



^ = b on dB 

dn 



for some /j, = (j,q > 0, where rj £ Cq°(B), < rj < 1 is a fixed radial cut-off function 
such that 

T](x) = 1 for \x\ < h, rj(x) = for \x\ > |. 

Step 2. Using a solution to (13) we construct, for some A > A*, a super-solution 
to (1). This provides a solution u\ for some A > A*, which is a contradiction. 

Proof of Step 1. Observe that given < R < 1 we must have for some Co = 
co (#) > 

(14) v(x) > u*(x) +c \x\<R. 

To prove this we recall the Green's function for A 2 with Dirichlct boundary condi- 
tions 

A 2 x G(x, y)=S y x £ B 
G(x, y)=0 x £ dB 
dC 

— (x,y) = x£dB, 

where S y is the Dirac mass at y £ B. Boggio gave an explicit formula for G(x,y) 
which was used in [16] to prove that in dimension N > 5 (the case 1 < N < 4 can 
be treated similarly) 

d{xfd{ V r 



(15) G(x,y)~\x-y\ 4 - N min[l 



\x - y\ 

where 

d(x) = dist(x,dB) = l-\x\ 



4 
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and a ~ b means that for some constant C > we have C a < b < Ca (uniformly 
for x, y € B). Formula (15) yields 

(16) G(x,y) >cd(x) 2 d(y) 2 

for some c > and this in turn implies that for smooth functions v and u such that 
i - u g ^(73) and A 2 (T) - u) > 0, 

ny) - m(w = y V 9n ^'^^ _ ~ A * & ( x >y) — ^ — J " x 

+ / G(x,y)A 2 (v - u) dx 
J B 

>cd(y) 2 [ (A 2 v - A 2 u)d{x) 2 dx. 
Jb 

Using a standard approximation procedure, we conclude that 

v(y) -u*(y) > cd{y) 2 X* [ (e v - e u * )d(x) 2 dx. 



Since v > u* , v ^ u* we deduce (14). 

Let u = (u* + v)/2. Then by Taylor's theorem 

(17) e° = e"° + (v - u )e u » + hv - u ) 2 e u ° + \{v - u ) 3 e u " + ±-(v - u fe^ 

I D 24 

for some Uq < £2 5= w and 
(18) 

e u * - e u ° + («* - u )e M ° + ^(u* - «o) 2 e" + «o) 3 e"° + ^(«* - ^o) 4 e 6 

2 o 24 

for some u* < £1 < Uo- Adding (17) and (18) yields 

(19) -{e v + e"*) > e M ° + -(« - u*) 2 e u °. 

2 8 

From (14) with R = 3/4 and (19) we see that uq = (u* + v)/2 is a super-solution 

of (13) with ^0 : = co/8. 

Proof of Step 2. Let us show now how to obtain a weak super-solution of (1) for 
some A > A*. Given /x > 0, let u denote the minimal solution to (13). Define tpi as 
the solution to 

A 2 <f 1 = firje u in B 
<Pi = on dB 

dip 



and if2 be the solution of 



„ - = on dB, 
an 



A 2 tp 2 = in B 

(f2 = a on dB 

dL P 2 u an 

— — = on dB. 

dn 



If JV > 5 (the case 1 < TV < 4 can be treated similarly), relation (16) yields 
(20) ipxix) > cxdix) 2 for all x e B, 
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for some c\ > 0. But u is a radial solution of (13) and therefore it is smooth in 
B \ -B1/4. Thus 

(21) u(x) < Mipt + ip 2 for all x £ B 1 / 2 , 

for some 71/ > 0. Therefore, from (20) and (21), for A > A* with A - A* sufficiently 
small we have 

(£ - l)u < <(>! + (£ - l)<f>2 inB. 

Let w = j^u — ^1 — (p- — 1)V2- The inequality just stated guarantees that w < u. 
Moreover 

k 2..._ x.u , V. 



and 



A 2 w = Ae" + — r\e u - me" > Ae" > Ae w in B 
A* 



w = a = b on dB. 

on 



Therefore w is a super-solution to (1) for A. By the method of sub and super- 
solutions a solution to (1) exists for some A > A*, which is a contradiction. □ 

Proof of Proposition 1.7. Let u £ H 2 (B), A > be a weak unbounded solution 
of (1). If A < A* from Lemma 2.6 we find that u < u\ where u\ is the minimal 
solution. This is impossible because u\ is smooth and u unbounded. If A = A* then 
necessarily u — u* by Theorem 1.2. □ 

Proof of Proposition 1.8. Let u denote the extremal solution of (1) with b > —4. 
If u is smooth, then the result is trivial. So we restrict to the case where u is 
singular. By Theorem 1.3 we have in particular that N > 13. We may also assume 
that a = 0. If b = -4 by Theorem 1.2 we know that if N > 13 then u = -4 log \x\ 
so that the desired conclusion holds. Henceforth we assume b > —4 in this section. 
For p > define 



so that 



Then 



u p (r) = u(pr) + 41og/9, 
A 2 u p = \*e u » in B 1/p . 



dp 

Hence, there is 5 > such that 



du n 

= u'(l)+4>0. 

p— l.r— 1 



u p (r) < u(r) for all 1 — S < r < 1, 1 — S < p < 1. 

This implies 

(22) u p (r)<u(r) for all < r < 1, 1 - 6 < p < 1. 
Otherwise set 

7*0 = sup{ < r < 1 I u p (r) > u(r) }. 

This definition yields 

(23) u p (r a ) = u(r ) and u' p (r )<u'(r ). 
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Write q = u(ro), j3 = u'(ro). Then u satisfies 





' A 2 u = 


= Ae 


(24) 


«(ro) = 


= a 




k "'M = 


= /?. 



(r) = u'(r)r + 4 for all < r < 1 

p=i 



on B ro 



Observe that u is an unbounded H 2 (B ro ) solution to (24), which is also stable. 
Thus Proposition 1.7 shows that u is the extremal solution to this problem. On the 
other hand u p is a supersolution to (24), since u' p (r ) < j3 by (23). We may now 
use Theorem 1.2 and we deduce that 

u(r) = u p (r) for all < r < r , 

which in turn implies by standard ODE theory that 

u(r) = u p {r) for all < r < 1, 

a contradiction with (22). This proves estimate (22). 
From (22) we see that 

rill 

(25) — (r) > for all < r < 1. 

Gt/9 p=l 

But 

and this together with (25) implies 

(26) ^r-(r) = -(u'{pr)pr + 4) > for all < r < -, < p < 1. 

dp p p 

which means that u p (r) is non-decreasing in p. We wish to show that lim p ^ u pi r ) 
exists for all < r < 1. For this we shall show 

(27) u p (r) > -41og(r) + log ^ 8 ( N ~ 2 )( N ~ for all < r < ^, < p < 1. 
Set 

'8 (AT -2)(N -Ay 
A* 

and suppose that (27) is not true for some < p < 1. Let 

ri = sup { < r < l/p | u p (r) < u (r) }. 

Observe that 

(28) A* > 8(7V — 2)(iV — 4). 

Otherwise w = — 41nr would be a strict supersolution of the equation satisfied by 
it, which is not possible by Theorem 1.2. In particular, r*i < l/p and 

u P {ri) = u (ri) and u' p (r\) > u'^n). 

It follows that mq is a supersolution of 



u o(r) = -41og(r) + log 



(29) 



' A 2 u 


= X*e 


< u 


= A 


du 






= B 


^ dn 



in B ri 
on dB ri 



on dB ri , 
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with A = Upiri) and B = u' (ri). Since u p is a singular stable solution of (29), it 
is the extremal solution of the problem by Proposition 1.7. By Theorem 1.2, there 
is no strict supcrsolution of (29) and we conclude that u p = uq first for < r < n 
and then for < r < 1/p. This is impossible for p > because u p (l/p) = 41ogp 
and u (l/p) < 41ogp + log( 8(jv ~ 2 A ) i jv ~ 4) ) < u p (l/p) by (28). This proves (27). 
By (26) and (27) we see that 

v(r) = lim u p (r) exists for all < r < +oo, 

p^O 

where the convergence is uniform (even in C k for any k) on compact sets of K Ar \{0}. 
Moreover v satisfies 

(30) A 2 v = X*e v in R N \ {0}. 
Then for any r > 

v(r) = lim u p (r) = lim u{pr) + 41og(pr) — 41og(r) = v(l) — 41og(r). 

p— »o p^O 

Hence, using equation (30) we obtain 

i f8(N -2)(N -4)\ 
w(r) = -41ogr + log( -i ^ '- J = « (r). 

But then 

u' p (r) = u'(pr)p — > —4, as p — > 0. 
and therefore, with r = 1 

(31) pu'(p)^-4 asp^O. 

□ 



4. Proof of Theorem 1.3 

We will show first 

Lemma 4.1. Suppose that the extremal solution u* to (1) is singular. Then for 
any a > there exists < R < 1 smc/i £/iai 

(32) (aO > (1 - cr) log f^), V \x\ < R. 

Proof. Assume by contradiction that (32) is false. Then there exists a > and a 
sequence x k £ B with x k — > such that 

(33) u*(a fc )<(l-<7)Iog^ 
Let Sfc = Ixfcl and choose < Afe < A* such that 



|x fc | 4 



(34) max u Xk = u Xk (0) = log f -j 

Note that \ k — > A*, otherwise m^,. would remain bounded. Let 

, , "A fc (sfcx) _ 1 

= "1 — 7TT" x € B k = —B. 
log(-i) Sfc 
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Then < v k < 1, v k (0) = 1, 

A 2 v k (x) = X.-^-e^^ 
log pr 

log(^-) 

by (34). By elliptic regularity v k — * w uniformly on compact sets of R w to a 
function v satisfying < v < 1, i>(0) = 1, A 2 v = in R w .By Liouville's theorem 
for biharmonic functions [17] we conclude that v is constant and therefore v = 1. 
Since \x k \ = s k we deduce that 

log(^) ' 

k 

which contradicts (33). □ 

Proof of Theorem 1.3. We write for simplicity u = u* , A = A*. Assume by 
contradiction that u* is unbounded and 5 < N < 12. If N < 4 the problem is 
subcritical, and the boundedncss of u* can be proved by other means : no singular 
solutions exist for positive A (sec [2])-though in dimension N = 4 they can blow up 
as A — > 0, see [24]. 

For e > let tp = \x\ ± ^ L+6 and let r} £ Cfi°(R N ) with r) = 1 in B 1/2 and 
supp(r]) C B. Observe that 

(A^) 2 = (F,v + 0(e))|x|- w+2e , where Hjy = 1 - V 

Using a standard approximation argument as in the proof of Lemma 2.6, we can 
use ipr) as a test function in (9) and we obtain 



(AV>) 2 + 0(1) > A / e u ^ 2 , 



B 

since the contribution of the integrals outside a fixed ball around the origin remains 
bounded as e — > (here O(l) denotes a bounded function as e — > 0). 
This implies 

(35) A / e-l^l 4 "™ < (H N + 0(e)) f \x\- N+2 ° = u> N ^ + 0(1) 

Jb Jb le 

where u>n is the surface area of the unit N — 1 dimensional sphere S 1 . In 
particular J B e u \x\ 4: ~ N+2e < +oo. 

For e > we define (f — \x\ 4 ~ N+2s . Note that away from the origin 

(36) A 2 ip = £k N \x\~ N+2e , where k N = A(N-2)(N -4) + 0(e). 
Let ifj solve 

f A 2 ^ = £fc Jv min(|x|- JV+2e , j) in B 

(37) ^ ^ ? 

ifj = -P- = on 9B. 



c9 
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Then ipj f ip as j — ► +00. Using (35) and (37) 

efcAT / u min(|x|^ Ar+2£ ,j) = / uA 2 ^- = A / e"^ 



< A / e> 



B 



<u N ^ + 0(l) 
2e 

where O(l) is bounded as e — > independently of j. Letting j —> +00 yields 

(38) ek N [ u \x\~ N+2e < uj n ^- + 0(1), 

Jb 2e 

showing that the integral on the left hand side is finite. On the other hand, by (32) 

(39) ek N f u \x\~ N+2e > ek N u N (l-a) f log(^)r- 1+2e dr = k N uj N (l - a)-. 

JB Jo r £ 

Combining (38) and (39) we obtain 

(l-a)k N <^ + 0(e). 



N 2 (N~4) 2 



Letting e — > and then a — > we have 

8(N-2)(N-i) <H N 

16 

This is valid only if N > 13, a contradiction. □ 

Remark 4.2. The conclusion of Theorem 1.3 can be obtained also from Proposi- 
tion 1.8. However that proposition depends crucially on the radial symmetry of the 
solutions, while the argument in this section can be generalized to other domains. 

5. The extremal solution is singular in large dimensions 

In this section we take a = b = and prove Theorem 1.4 for N > 32. 

The idea for the proof of Theorem 1.4 is to to estimate accurately from above the 
function A*e M , and to deduce that the operator A 2 — A*e" has a strictly positive 
first eigenvalue (in the Hq(B) sense). Then, necessarily, u* is singular. 

Upper bounds for both A* and u* are obtained by finding suitable sub and 
supersolutions. For example, if for some Ai there exists a supersolution then A* > 
Ai. If for some A2 one can exhibit a stable singular subsolution u, then A* < A2. 
Otherwise A2 < A* and one can then prove that the minimal solution u\ 2 is above 
u, which is impossible. The bound for u* also requires a stable singular subsolution. 

It turns out that in dimension N > 32 we can construct the necessary subso- 
lutions and verify their stability by hand. For dimensions 13 < N < 31 it seems 
difficult to find these subsolutions explicitly We adopt then an approach that 
involves a computer assisted construction of subsolutions and verification of the 
desired inequalities. We present this part in the next section. 

Lemma 5.1. Assume N > 13. Then u* < u = — 41og|x| in B\. 
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(40) 



Proof. Define u(x) = — 41og|x|. Then u satisfies 

A 2 u = 8(JV-2)(JV-4)e° in R N 
u = on dB\ 

^ = -4 on 8B 1 

on 

Observe that since u is a supersolution to (1) with a = b = we deduce imme- 
diately that A* > 8(JV - 2)(JV - A). 

In the case A* = 8(JV — 2){N — 4) we have u\ < u for all < A < A* because u is 
a supersolution, and therefore u* < u holds. Alternatively, one can invoke Theorem 
3 in [2] to conclude that we always have A* > 8(JV - 2)(JV — 4). 

Suppose now that A* > 8(JV - 2){N - 4). We prove that u x < u for all 8(JV - 
2) (JV — 4) < A < A* . Fix such A and assume by contradiction that u\ < u is not 
true. Note that for r < 1 and sufficiently close to 1 we have u\(r) < u(r) because 
u' x {l) = while u'{\) = -4. Let 

Ri = inf{ < R < 1 | u x < u in (R, 1) }. 

Then < R\ < 1, u\(R\) = u(Ri) and u' x (Ri) < u'(Ri). So u\ is a super-solution 
to the problem 

A 2 u = 8(N - 2) (JV - 4)e" in 

u = u\(R 1 ) ondB Rl 
dii 

— =u' x {Ri) ondB Rl 

while u is a subsolution to (40). Moreover it is stable for this problem, since from 
Rellich's inequality (6) and 8(JV-2)(JV-4) < N 2 (N - 4) 2 /16 for JV > 13, we have 

8(JV-2)(JV-4) / eV< ^ (j !~ 4)3 / t4</ (A^) 2 V^C °°(%). 

By Remark 2.7 we deduce u < u\ in Bjjj which is impossible. □ 

An upper bound for A* is obtained by considering again a stable, singular sub- 
solution to the problem (with another parameter, though): 

Lemma 5.2. For N > 32 we have 

(41) A* < 8(JV - 2)(JV - 4)e 2 . 

Proof. Consider w = 2(1 — r 2 ) and define 



where u(x) = — 41og|x|. Then 

A 2 u = 8(JV - 2)(JV - 4)4 = 8(JV - 2)(JV - 4)e a = 8(JV - 2)(JV - A)e u+W 
r 4 

<8(JV-2)(JV-4)e 2 e". 

Also u(l) = u'(l) — 0, so u is a subsolution to (1) with parameter Ao = 8(JV — 
2)(JV-4)e 2 . 

For JV > 32 we have A < JV 2 (JV - 4) 2 /16. Then by (6) u is a stable subsolution 
of (1) with A = A . If A* > A = 8(JV-2)(JV-4)e 2 the minimal solution u X[> to (1) 
with parameter Aq exists and is smooth. From Lemma 2.6 we find u < u\ which 
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is impossible because u is singular and u\ is bounded. Thus we have proved (41) 
for N > 32. □ 

Proof of Theorem 1.4 in the case N > 32. 

Combining Lemma 5.1 and 5.2 we have that if N > 32 then X*e u ' < r~ 4 8(A - 
2)(N - 4)e 2 < r- i N 2 (N - 4) 2 /16. This and (6) show that 

which is not possible if u* is bounded. □ 



6. A COMPUTER ASSISTED PROOF FOR DIMENSIONS 13 < N < 31 



Throughout this section we assume a = b = 0. As was mentioned in the previous 
section, the proof of Theorem 1.4 relies on precise estimates for u* and A*. We 
present first some conditions under which it is possible to find these estimates. 
Later we show how to meet such conditions with a computer assisted verification. 

The first lemma is analogous to Lemma 5.2. 

Lemma 6.1. Suppose there exist e > 0, A > and a radial function u € H 2 {B) n 



„ (B \ {0}) such that 

A 2 u < Xe" 



\u(l)\<s, 
u # L°°{B) 



for allO < r < 1 
du , 



On 



(1) 



< e 



(42) 
Then 

Proof. Let 

(43) 
so that 

and 

It follows that 



Ac 



e if 



< 



(Af) 2 for all(f€C%°(B). 



A* < Xe 



2s 



ip(r) 



er 



2s 



A 2 V> = 0, ip(l) = -e, V'(l) = 2e 
-2s < ip(r) < -e for all < r < 1. 
A 2 (u + ip)< Xe u = Xe^e^ < Xe 2e e u+ ^. 



On the boundary we have u(l) + "0(1) < 0, u'(l) + ip'(l) > 0. Thus u + ip is a 
singular subsolution to the equation with parameter Ae 2£ . Moreover, since ip < — £ 
we have Xe 2e e u+ ^ < Xe £ e u and hence, from (42) we see that u + ip is stable for the 
problem with parameter Ae 2e . If Ae 2e < A* then the minimal solution associated to 
the parameter Ae 2e would be above u + ip, which is impossible because u is singular. 
□ 
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Lemma 6.2. Suppose we can find e > 0, A > and u e H 2 (B) n W£™(B \ {0}) 
such that 

A 2 u > Xe u for all < r < 1 



l«(l)l<e, 



ITU 

on 



< e. 



Then 



Ac 



-2s 



< A* 



Proof. Let ip be given by (43). Then u — i\) is a supersolution to the problem with 
parameter Xe~ 2e . □ 

The next result is the main tool to guarantee that u* is singular. The proof, as 
in Lemma 5.1, is based on an upper estimate of u* by a stable singular subsolution. 

Lemma 6.3. Suppose there exist eo,s > ; A a > and a radial function u € 
H 2 (B) n Wf'°°(B \ {0}) such that 



(44) 


A 2 u < (A a 


+ 


eo)e u 


for all < 


r < 1 


(45) 


A 2 u > (A a 




eo)e a 


for all < 


r < 1 


(46) 


Kl)l <e, 




> 


< e 




(47) 


u L°°(B) 










(48) 


/9b / eV 
Jb 


< 


[ (Atp) 2 for all 
Jb 


<P G C °°(S), 


where 












(49) 




a (A a + S ) 9e 
(A a - Sq) 2 




Then u* 


is singular and 










(50) 


(A„- 




)e~ 2e < X* < (A Q + 


£o)e 2£ . 



Proof. By Lemmas 6.1 and 6.2 we have (50). Let 

A a + £0 



f) = log 



A a — £0 

S A 



lie. 



and define 

ip(r) =-"-r* + 28. 
We claim that 

(51) u* < u + ip in Bi. 
To prove this, we shall show that for A < A* 

(52) u\ < u + tp in B\. 
Indeed, we have 

A 2 ip = -S2N(N + 2) 

ip(r) > S for all < r < 1 

(p(l)>5>e, <p'(l) = -5< 
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and therefore 

A 2 (u + <p)< (A„ + e )e u + A 2 ^ < (A. + e Q )e u = (X a + e^e^ 

(53) < (X a +e Q )e- s e u+v . 
By (50) and the choice of 5 

(54) (A„ + e )e- 5 = (A a - £oK 3£ < A*. 

To prove (52) it suffices to consider A in the interval (A a — £o)e -3e < A < A*. Fix 
such A and assume that (52) is not true. Write 

u = u + tp 

and let 

Ri = sup{ < R < 1 1 u x (R) = u(R) }. 
Then < R\ < 1 and u x {Ri) = u{Ri). Since u' x (l) = and u'(l) < we must 
have u' x (Ri) < u'(Ri). Then u\ is a solution to the problem 

A 2 u = \e u in J B 7?1 
u = u\(Ri) on dB Rl 
3u 

— = u x (Ri) on dB Rl 
an 

while, thanks to (53) and (54), u is a subsolution to the same problem. Moreover 
u is stable thanks to (48) since, by Lemma 6.1, 

(55) A < A* < (A + e )e 2e 
and hence 

Ae a < (Xa+eo)e 2£ e 2S e u <p e u . 

We deduce u < u\ in B Rl which is impossible, since u is singular while u\ is 
smooth. This establishes (51). 
From (51) and (55) we have 

A*e u * < /3 e- E e u 

and therefore 

f'(Aip) 2 — X*e u * ip 2 
inf JbK , — > 0. 

^Cg-(B) J B p 2 

This is not possible if u* is a smooth solution. □ 

For each dimension 13 < N < 31 wc construct u satisfying (44) to (48) of the 
form 



(56) u(r) = 



-41ogr + log( 8 ^- 2 ^- 4) for0<r<r 



A 

^u(r) for tq < r < 1 

where w is explicitly given. Thus u satisfies (47) automatically. 
Numerically it is better to work with the change of variables 

w(s) = u(e s ) + 4s, -oo < s < 

which transforms the equation A it = Xe u into 

Lw + 8(N -2)(N - 4) = Xe w , -oo < s < 
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where 

£ + ^ - ^ + ^ - ™ + 2 4? - ^ -^-^ 

The boundary conditions u(l) = 0, u'(l) = then yield 

w(0)=0, w'(0)=4. 

Regarding the behavior of wass-t — oo observe that 

, /8(iV-2)(iV-4)\ , 
u[r) = — 41ogr + log I I for r < r 



if and only if 



, 8(iV-2)(iV-4) 

ii>(s) = log for all s < logr . 

A 



The steps we perform are the following: 

1) We fix xo < and using numerical software we follow a branch of solutions to 

' Lw + 8{N -2) (AT -4) = \e™, x < s < 

w(0) = 0, w'(0)=t 

8(AT-2)(iV-4) d 2 w d 3 w 
w{x )=log , —^-{x ) = 0, -—p(xo) = 

as t increases from to 4. The numerical solution (w, A) we are interested in 
corresponds to the case t = 4. The five boundary conditions are due to the fact 
that we are solving a fourth order equation with an unknown parameter A. 

2) Based on w, A we construct a C 3 function w which is constant for s < Xq 
and piecewise polynomial for xq < s < 0. More precisely, we first divide the 
interval [xo,0] in smaller intervals of length h. Then we generate a cubic spline 
approximation gji with floating point coefficients of 4-!r- From gji we generate a 
piecewise cubic polynomial g ra which uses rational coefficients and we integrate it 
4 times to obtain w, where the constants of integration are such that ^t(x ) = 0, 
f < j < 3 and w(x n ) is a rational approximation of log(8(7V — 2) (N — 4) /A). Thus w 
is a piecewise polynomial function that in each interval is of degree 7 with rational 
coefficients, and which is globally C 3 . We also let A be a rational approximation of 
A. With these choices note that Lw + 8(N — 2)(N — 4) — Xe w is a small constant 
(not necessarily zero) for s < xq. 

3) The condition (44) and (45) we need to check for u are equivalent to the 
following inequalities for w 

(57) Lw + 8(N - 2){N - 4) - (X + e a )e w < 0, -oo < s < 

(58) Lw + 8(N -2)(N -4)- (\-e )e w > 0, -oo < s < 0. 

Using a program in Maple we verify that w satisfies (57) and (58). This is done 
evaluating a second order Taylor approximation of Lw + 8(N — 2)(N — 4) — (X + e)e w 
at sufficiently close mesh points. All arithmetic computations are done with rational 
numbers, thus obtaining exact results. The exponential function is approximated by 
a Taylor polynomial of degree 14 and the difference with the real value is controlled. 
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More precisely, we write 

f(s) = Lw + 8(N ~ 2)(N - 4) - (A + e )e w , 

f(s) =Lw + 8(N - 2)(N - 4) - (A + e )T(w), 

where T is the Taylor polynomial of order 14 of the exponential function around 0. 
Applying Taylor's formula to / at yj, we have for s € [yj, yj+h], 

f{s)<f{y 3 ) + \.f'{y 3 )\h+ l -Mh 2 



< f(3/i) + \f'(yj)\h + \Mh 2 + \f( Vj ) - f( yj )\ + \f'(y 3 ) - f'( yj )\h 

< Kvi) + \f'(y 3 )\ h + l Mh2 + E i + E 2K 



where 



M is a bound for |/"| in [yj,yj + h] 

E x is such that (A + e )\e w - T(w)\ < E x in [yj,yj + h] 

E 2 is such that (A + e )\{e w - T'(w))w'\ < E 2 in [yj,yj + h}. 

So, inequality (57) will be verified on each interval [yj, yj + h] where w is a polyno- 
mial as soon as 

(59) f{Vj) + \f\Vj)\h + \Mh 2 + E,+ E 2 h < 0. 

When more accuracy is desired, instead of (59) one can verify that 

Rxi) + \f'(xi)\- + hl(-) 2 + E 1+ E 2 -<0, 
m 2 m m 

where (xi)i—i,,, m +i are m + 1 equally spaced points in [yj, yj + h}. 

We obtain exact values for the upper bounds M, E\, E 2 as follows. First note that 

/" = W-(A + e )e^((u;') 2 + ^")- On [y 3 ,y 3 +h], we have w(s) = ELo^( s -%) 4 

and we estimate |w(s)| < Y^l=o I *!' 1 * f° r s G \Vii Vj + h]. Similarly, 



(60) 



At 7 

aw 



<^2i(i-l)...(i-£ + 1)1^:1^ for all s e [yj, Vj + h] 



The exponential is estimated by e w < e 1 < 3, since our numerical data satisfies the 
rough bounds —3/2 < w < 1. Using this information and (60) yields a rational 
upper bound M. E\ is estimated using Taylor's formula : 

(3/2) 15 



E 1 = (A + £ )- 



15! 



Similarily, E 2 = (A + £0) ^(41 B\ where B\ is the right hand side of (60) when 



14! 
1=1. 



4) We show that the operator A 2 — (3e u where u(r) = w(logr) — 41ogr, satisfies 
condition (48) for some (3 > (3q where f3o is given by (49). In dimension N > 13 
the operator A 2 — f3e u has indeed a positive eigenfunction in H 2 (B) with finite 
eigenvalue if (3 is not too large. The reason is that near the origin 
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where c is a number close to 8(JV — 2)(JV— 4)/3/A. If (3 is not too large compared to A 
then c < N 2 (N — 4) 2 /16 and hence, using (6), A 2 — j3e u is coercive in HQ(B ro ) (this 
holds under even weaker conditions, see [11]). It follows that there exists a first 
eigcnfunction (pi e Hq(B) for the operator A 2 — (3e u with a finite first eigenvalue 
/ii, that is 

A 2 (^i — f3e u ipi = [iipi in B 
tpi > in B 
<Pi G H*{B). 
Moreover fi\ can be characterized as 



/i! = inf JB 



/ B (A^) 2 -/3e> 



2 



<pecg°(B) f B ip 2 

and is the smallest number for which a positive eigcnfunction in Hq(£1) exists. 

Thus to prove that (48) holds it suffices to verify that [i\ > and for this it is 
enough to show the existence of a nonnegative tp <G H 2 (B), <p ^ such that 

( A 2 ip-f3e u <p > inB 
tp = on dB 



(61) 



^ < on dB. 
on 



Indeed, multiplication of (61) by ipi and integration by parts yields 

Mi / f<Pi + ( ^f~Aip 1 > 0. 
Jb JdB o n 

But Aipi > on dB and thus /ii > 0. To achieve (61) we again change variables 
and define 

<fi(s) = ip(e s ) - oo < s < 0. 
Then we have to find <j) > 0, <j) ^ satisfying 

' £c/> - ^e 1 "^ > in -oo < s < 
(62) J 0(0) = 

«//(0) < 0. 

Regarding the behavior as s — » — oo, we note that i« is constant for — oo < s < xq, 
and therefore, if 

L<f> - /Je 1 "^ = - oo < s < x 

then <f> is a linear combination of exponential functions e~ as where a must be a 
solution to 

a 4 - 2(JV - 4)a 3 + (TV 2 - 107V + 20)a 2 + 2(JV - 2)(JV - 4) a = /Je™^ 

where ^e^* ) is close to 8(JV - 2)(JV - 4)/3/A. If JV > 13 the polynomial 

a 4 - 2{N - 4)a 3 + (JV 2 - 107V + 20)a 2 + 2(JV -2)(N- 4)q - 8(JV - 2)(JV - 4) 

has 4 distinct real roots, while if JV < 12 there are 2 real roots and 2 complex 
conjugate. If JV > 13 there is exactly one root in the interval (0, (JV — 4)/2), 2 roots 
greater than (JV — 4)/2 and one negative. We know that ip(r) = </>(logr) ~ r~ a is 
in H 2 , which forces a < (JV — 4)/2. It follows that for s < xq, <j> is a combination of 
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e~ a ° s , e _QlS where ao > 0, a\ < are the 2 roots smaller than a < (N — 4)/2. For 
simplicity, however we will look for </> such that <fi(s) = Ce~ a ° s for s < xo where 
C > is a constant. This restriction will mean that we will not be able to impose 
4>'(0) — at the end. This is not a problem because <fi'(0) < 0. 

Notice that we only need the inequality in (62) and hence we need to choose 
a G (0, N - 4/2) such that 

a 4 - 2(N - 4)a 3 + (N 2 - 107V + 20)a 2 + 2(N - 2)(N - 4)a > (3e w{xo) . 

The precise choice we employed in each dimension is in a summary table at the end 
of this section. 

To find a suitable function (f> with the behavior (f>(s) = C'e~ a for s < xo we set 
4> = Tpe~ as and solve the following equation 

T a ^ - (3e w ^j = f 

where the operator T a is given by 

^4 / j3 / j2 

Ta ^ = W + ^ a + 2(7V ~ 1£ + ^ ~ 6a(7V _ ^ + ~ WN + 20) ^2" 

+ (-4a 3 + 6a 2 (N - 4) - 2a(N 2 - ION + 20) - 2(N - 2)(N - 4))-^ 

ds 

+ (a 4 - 2a 3 (N - 4) + a 2 (N 2 - 10 AT + 20) + 2a(N -2)(N- 4))^ 

and / is some smooth function such that / > 0, / ^ 0. Actually we choose (3 > /3q 
(where /3o is given in (49)) find a satisfying approximately 

a 4 - 2(N - 4)a 3 + (N 2 - 10N + 20)a 2 + 2 (AT - 2) (AT - 4)a = (3e w{xo) . 

We solve numerically 

- (3e w i> = f x <s<0 
$(x ) = 1, $"(xo) = 0, 4>"'(x ) = 
^(0) = 

Using the same strategy as in 2) from the numerical approximation of -j-r we 
compute a piecewise polynomial ip of degree 7, which is globally C 3 and constant 
for s < xo. The constant ip(xo) is chosen so that ip(0) = 0. We use then Maple to 
verify the following inequalities 

ip > x < s < 
T Q V - /Se^V > x < s < 
V>'(0) < 

where (3o < (3 < (3 and < a < (N — 4)/2 are suitably chosen. 

At the URLs: 
http : //www . lamf a . u-picardie . f r/ dupaigne/ 
http : //www. ime .unicamp.br/~msm/ 

we have provided the data of the functions w and tp defined as piecewise polynomials 
of degree 7 in [a;o,0] with rational coefficients for each dimension in 13 < A" < 31. 
Wc also give a rational approximation of the constants involved in the corresponding 
problems. 
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We use Maple to verify that w and ip (with suitable extensions) are C 3 global 
functions and satisfy the corresponding inequalities, using only its capability to 
operate on arbitrary rational numbers. These operations are exact and are limited 
only by the memory of the computer and clearly slower than floating point opera- 
tions. We chose Maple since it is a widely used software, but the reader can check 
the validity of our results with any other software (see e.g. the open-source solution 
pari/gp). 

The tests were conducted using Maple 9. 
Summary of parameters and results 



N 


A 


£o 




£ 







P 


a 


13 


2438.6 


1 


5 


10" 


-7 


2550 


2500 


3.9 


14 


2911.2 


1 


3 


10 


-6 


3100 


3000 


3.4 


15 


3423.8 


1 


3 


10 


-6 


3600 


3500 


3.1 


16 


3976.4 


1 


1 


10 


-5 


4100 


4000 


3.0 


17 


4568.8 


1 


2 


10 


-4 


4800 


4600 


3.0 


18 


5201.1 


2 


2 


10 


-4 


5400 


5300 


2.7 


19 


5873.2 


2 


2 


10 


-4 


6100 


6000 


2.7 


20 


6585.1 


3 


7 


10 


-4 


7000 


6800 


2.7 


21 


7336.7 


3 


7 


10 


-4 


7700 


7500 


2.6 


22 


8128.1 


4 


1 


10 


-3 


8600 


8400 


2.6 


23 


8959.1 


4 


1 


10 


-3 


9400 


9200 


2.5 


24 


9829.8 


4 


1 


10 


-3 


10400 


10200 


2.5 


25 


10740.1 


4 


1 


10 


-3 


11400 


11200 


2.5 


2G 


11690.1 


6 


2 


10 


-3 


12400 


12200 


2.5 


27 


12679.7 


7 


2 


10 


-3 


13400 


13200 


2.4 


28 


13709.0 


7 


2 


10 


-3 


14500 


14300 


2.4 


29 


14777.8 


7 


2 


10 


-3 


15400 


15200 


2.4 


30 


15886.2 


8 


2 


10 


-3 


16600 


16400 


2.4 


31 


17034.3 


10 


2 


10 


-3 


17600 


17500 


2.3 



Remark 6.4. 1 ) Although we work with A rational, in the table above we prefer to 
display a decimal approximation of X. 

2) In the previous table we selected a "large" value of eo in order to have a fast 
verification with Maple. By requiring more accuracy in the numerical calculations, 
using a smaller value of £o and using more subintervals to verify the inequalities 
in the Maple program it is possible to obtain better estimates of A* . For instance, 
using formulas (50), we obtained 

N ^ £0 £ Knm Knax (3 (3 a 

13 2438.589 0.003 5-10" Y 2438.583 2438.595 2550 2510 3.9 

14 2911.194 0.003 5-10" 7 2911.188 2911.200 3100 3000 3.4 

The verification above, however, required to check 1500 subintervals of each of 
the 4500 intervals of length 0.002, which amounts to substantial computer time. 

7. Proof of Proposition 1.5 

Throughout this section, we restrict, as permitted, to the case a = 0. 
a) Let u denote the extremal solution of (1) with homogeneous Dirichlet boundary 
condition a = b = 0. We extend u on its maximal interval of existence (0, R). 
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Lemma 7.1. R < oo and u(r) ~ log(i? — r) 4 /or r ~ i?. 

Proof. The fact that R < oo can be readily deduced from Section 2 of [1]. We 
present an alternative (and more quantitative) argument. We first observe that 



(63) 



i"--u'>0 for all r G [1, R). 



Integrate indeed (1) over a ball of radius r to conclude that 



(64) 



< A 



d A 
— Au = w N r 

or 



N-l I u m 



N-l , „ 1 
(u u ) 

r r 



If r = 1, since u is nonncgativc in (0,1) and u(l) = = 0, we must have 

u"(l) > 0. In fact, u"(l) > 0. Otherwise, we would have u"{l) = and u"'(l) > 
by (64), contradicting u > in (0, 1). So, we may define 

R = sup{r > 1 : u" (t) - -u'(t) > for all t G [1, r)} 

and we just need to prove that R = R. Assume this is not the case, then u"(R) 



(u» - £«' 



)' (R) < 0. This contradicts (64) and we have 



^u'(R) = and u"'(R) 
just proved (63). In particular, we see that u is convex increasing on (1, R). 



Since u is radial, (1) reduces to 
„(4) + 2(7V-1) 



(65) 
Multiply by it' 



(JV-l)(JV-3) . 7< _ (JV-l)(JV-3) .., _ 

,-, Lt „ IX /AC- 



u(V + 2 ^- 1 ) u »v + 

r 



, (N -l)(N-3) „ , (N-l)(N-3) , , a , 
5 u w ; (u ) = A(e ) , 



which we rewrite as 

[(u'"u'Y - u"'u"} + 2{N - 1) 



+ {N-l)(N-3) 



2r 2 



A(e")' 



By (63), it follows that for r £ [1, R), 

1 



[{u"'u')' - u"'u"\ + 2{N - 1) -u"u' + (N - l)(N - 3) 



2r 2 



> A(e u )' 



Integrating, we obtain for some constant A 



(«") 



1 



+ 2(7V- l)-u"u' 
2 r 



(JV-l)(JV-3) {u'f 



> Ae" - A. 



We multiply again by v! 
(66) 



(u'V) 2 )' ~ «" (K) 2 )'l - -{u"fu' + 2(iV - l)-u"(u') 2 
x ' \ 2 r 



We deduce from (63) that 



1 



u"{u'f 



1 / 1 



2 



'\3 



:(«) 



( W ') 2 (-')' < J ( iK) 3 
r 2 \ r 



and 
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V) 3 < V) v' < i ( V) 3 

r z r 2 \r 

Using this information in (66), dropping nonpositive terms and integrating, we 
obtain for some constant B, 

u"(u') 2 + ^—^-(u') 3 > Xe u -Au-B 
4 r 

Applying (63) again, it follows that for C = 4 -1 + 1 

Cu"{u') 2 >Xe u -Au-B 

which after multiplication by u' and integration provides positive constants c, C 
such that 

{u'f > c(e u - Au 2 -Bu- C). 

At this point, we observe that since u is convex and increasing, u converges to +oo 
as r approaches R. Hence, for r close enough to R and for c > perhaps smaller, 

u' > c e u/i . 

By Gronwall's lemma, R is finite and 

u < —4 \og(R — r) + Cfor rclose to R. 

It remains to prove that u > — 41og(-R — r) — C. This time, we rewrite (1) as 

y-^Au)']' = \r Jf - 1 e u . 
We multiply by r JV_1 (Au)' and obtain : 

X - [r 2N - 2 {{Au)') 2 }' = Xr 2N - 2 e u {Au)' < Ce u (Au)' < C{e u Au)'. 

Hence, for r close to R and C perhaps larger, 

((Am)') 2 < Ce u Au 

and so 

VAH(Au)' < Ce u/2 Au < C'e u/2 u" < C'(e u/2 u'Y, 
where we have used (63). Integrate to conclude that 

(Au) 3 / 2 < Ce"/V. 
Solving for Au and multiplying by (u') 1 ' 3 , we obtain in particular that 

(u') 1/3 u" < Ce"/V. 
Integrating again, it follows that (u') 4 / 3 < Ce u ^ 3 , i.e. 

u' < Ce u/A . 
It then follows easily that (for r close to R) 

u > -Alog(R-r) - C. 

□ 

Proof of Proposition 1.5 a). Given N > 13, let b max denote the supremum of 
all parameters b > —4 such that the corresponding extremal solution is singular. 
We first observe that 

b raax > ^ 
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In fact, it follows from Sections 5 and 6 that the extremal solution u associated to 
parameters a = b = is strictly stable : 

67 inf jBy ^' — — > 0. 

Extend u as before on its maximal interval of existence (0, R). Choosing R £ (1, R) 
close to 1, we deduce that (67) still holds on the ball Br. In particular, letting 
v(x) — u(Rx) — u(R) for x 6 B, we conclude that v is a singular stable solution 
of (1) with a = and b = Ru'(R) > 0. By Proposition 1.7, we conclude that 
b max > 0. We now prove that 

b max < oo. 

Assume this is not the case and let u n denote the (singular) extremal solution 
associated to b n , where b„ / oo. We first observe that there exists p n £ (0, 1) such 
that u' n {p n ) = 0. Otherwise, u n would remain monotone increasing on (0, 1), hence 
bounded above by u n (l) — 0. It would then follow from (1) and elliptic regularity 
that u n is bounded. Let v n (x) = u n {p n x) — u n {p n ) for x € B and observe that v n 
solves (1) with a = b = and some A = A„. Clearly v n is stable and singular. By 
Proposition 1.7, v n coincides with u, the extremal solution of (1) with a = b = 0. 
By standard ODE theory, v n — u on (0, R). In addition, 



1 

Pn 





= 1.' 




(p n ) 




On) 




Pn 





-co, 



which can only happen if l/p n — > R- 

Now, since u n is stable on B, u — v n is stable on B\/ Pn . Letting n — * oo, we 
conclude that u is stable on B^. This clearly contradicts Lemma 7.1. 

We have just proved that b max is finite. It remains to prove that u* is singular 
when —A < b < b max . We begin with the case b = b max . Choose a sequence 
(b n ) converging to b max and such that the corresponding extremal solution u n is 
singular. Using the same notation as above, we find a sequence p n € (0, 1) such 
that 

1 ,/ 1 
— u 



Pn \Pn, 

Taking subsequences if necessary and passing to the limit as n — » oo, we obtain for 
some p € (0, 1) 



1 



- =b 



III (IX 



— U 

P \P; 

Furthermore, by construction of p n , u is stable in B\/ Pn hence in B 1 / p . This implies 
that v defined for x £ B by v(x) = u(-) — u(-), is a stable singular solution of 
(1) with b = b max . By Proposition 1.7, we conclude that the extremal solution is 
singular when b = b max . 

When b = —4 we have already mentioned in the introduction that u* is singular 
for N > 13 as a direct consequence of Proposition 1.7 and Rellich's inequality. 

So we are left with the case —4 < b < b max . Let u* m denote the extremal 
solution when b = b max , which is singular, and A^ the corresponding parameter. 
For < R < 1 set 

u R (x) = u* m (Rx) - u* m (R). 
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Then 

A 2 u R = X R e UR where X R = A£i? 4 e"™ (ii) , 
and ur = on dB, while 

du R _ u du* 

By (31), note that 

du* du* 
R^(R)^b max aaR->l, and R—^(R) -> -4 as R -> 0. 
ar ar 

Thus, for any —4 < b < b max we have found a singular stable solution to (1) (with 
a = 0). By Proposition 1.7 the extremal solution to this problem is singular. 

□ 

Proof of Proposition 1.5 b). Let b > —4. Lemma 5.1 applies also for b > —4 
and yields u* < u where u(x) = — 41og|ar|. We now modify slightly the proof of 
Lemma 5.2. Indeed, consider w = (4 + b)(l — r 2 )/2 and define u = u — w. Then 

A 2 u = 8(N - 2)(N - 4)\ = 8(N - 2)(N - 4)e" = 8(N - 2)(N - A)e ll+W 

< 8(7V -2){N -A)e^ +b ^ 2 e u . 

Also u(l) = 0, u'(l) = b, so u is a subsolution to (1) with parameter Ao = 8(N — 
2)(7V-4)e( 4+b )/ 2 . 

If N is sufficiently large, depending on b, we have Ao < N 2 (N — 4) 2 /16. Then 
by (6) u is a stable subsolution of (1) with A = Ao- As in Lemma 5.2 this implies 
A* < A . 

Thus for large enough TV we have A*e u * < r" 4 8(iV - 2) (TV - 4)e( 4+b '/ 2 < 
r~ 4 N 2 (N - 4) 2 /16. This and (6) show that 

v>ec~(B) \ B ip 2 
which is not possible if u* is bounded. □ 
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